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We study the dynamics of two particles that interact only when in contact. In this sense, although 
not in every particular, the interactions mimic those in granular materials. The detailed solution of 
the dynamics allows an analysis of the backscattering behavior of the first particle and of the energy 
dissipation in the system as a function of various parameters. 
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I. INTRODUCTION 

Signal propagation in granular materials is a long-standing subject of interest 0,0,111, 0,0], an d has recently attracted 
considerable attention in a number of novel contexts. One is the observation that the attenuation characteristics of 
signals (impacts) in such materials might make them good candidates for shock absorpti on Ifj . Another is the use of 
the backscattered signal for the nondestructive identification of buried objects 0, 0> B 0> El- 
Interactions in granular materials are notoriously complex, but a number of relatively simple models have been 
implemented in the context of signal propagation, models that it is hoped capture the principal ingredients of the 
true interactions at least in some regimes. A model potential used for a monodisperse chain of particles that repel 
upon overlap according to the Hertz law is given by 

V(S iA+1 ) =aSf ti+1) S<0, 

(1) 

V(Si, i+1 ) =0, 5>0. 

Here 

5i,i+i =2R - [{zi+i + x i+ i) - fa + Xi)] , (2) 



Y and a denote Young's modulus and Poisson's ratio, Zi denotes the initial equilibrium position of grain i in the 
chain, and xt is the displacement of grain i from this equilibrium position. The geometric parameter R is the radius if 
the particles are elastic spheres. More generally, R is determined by the principal radii of curvature of the surfaces at 
the point of contact |12|. The exponent n is 5/2 for spheres, it is 2 for cylinders, and in general depends on geometry. 

Two features of the potential are decisive for the resulting propagation properties of the chain: one is the exponent 
n, and the other is the absence of an attractive portion in the potential. It is in this latter respect that the granular 
chain differs profoundly from the many models of nonlincarly interacting oscillators (such as the Fermi-Pasta-Ulam 
chain) extensively studied in recent years [l3L Il4| . 

Two decades ago Nesterenko showed that under appropriate assumptions, among them the slow spatial variation 
of the displacements xt , the equations of motion for the granular particles could be approximated by a continuous 
nonlinear partial differential equation that admits a soliton solution for a propagating perturbation in the chain. 
More recently it was shown that strictly speaking these solutions are solitary waves rather than solitons [Hill. It 
is also understood that a potential exponent n > 2 is required to accommodate solitary waves, i.e. that a "hard 
nonlinearity" is required Q- Nesterenko coined the words "sonic vacuum" to stress the importance of the absence 
of "harmonic terms" (n — 2) in this analysis because such terms (and only such terms) are associated with ordinary 
sound propagation. Indeed, as n — > 2 the width of the solitary wave diverges. More recently, numerical simulations of 
discrete chains have provided a great deal of additional information about the important role of discreteness as well 
as nonlinearity in the signal propagation properties. 

Thus, while the continuum solitary wave solutions provide enormously helpful analytic insights, the importance of 
discreteness in the propagation properties of such chains should not be underestimated. Neither should the fact that 
even with a "harmonic" exponent n — 2, the chain is decidedly different from a harmonic system because there is no 
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attractive restoring force. The resulting nonanalyticity of the potential leads to profoundly different properties than 
those observed in an ordinary harmonic chain. Indeed, although solitary waves are no longer exact solutions (and 
there is no longer a perfect sonic vacuum), there are extremely long lived quasi- solitary waves that broaden slowly 
(whereas a harmonic chain would exhibit rapid dispersion of energy) [l4j |. 

We offer the above observations in support of the usefulness of a study of a chain with only repulsive forces even 
with the choice n — 2. A number of helpful insights can be achieved from such a study if one is able to obtain analytic 
results, and this is at all possible only with n = 2. In this paper we carry out such an analysis for the simplest 
(and yet extremely interesting) such system, namely, one consisting of only two masses. Surprisingly, this simple 
scenario seems not to have been analyzed in detail in the literature, although the dynamics of pairwise collisions is 
an important input (and phenomenological assumptions are made about it) in the study of longer chains ||. A two- 
mass system, being explicitly integrable, can yield valuable information. In view of some of the current applications 
of energy propagation in granular materials, we are particularly interested in backs cattering and energy dissipation 
in the presence of friction, and on the effects of frictional disorder on these quantities. The results of our analysis 
promise to provide helpful phenomenological input into the propagation of excitation pulses in longer chains as long 
as the pulses remain narrow. For example, recent experiments on impulse dispersion in tapered chains indicate the 
significant role of friction in the propagation and dispersion process |l7j . We find that a two-mass system leads to 
some unanticipated nonmonotonic (albeit, in hindsight understandable) behaviors as one varies system parameters. 
To support the relevance of our results to the spherical granular problem, we also present some numerical results for 
the Hertz potential with n = 5/2 and show the similarity in the behavior of this system and that of the n — 2 case. 

We consider two finite-sized particles ("granules") of equal mass m. The specific value of the mass does not matter 
since it can be scaled out by redefining variables; the only important point is that the masses are equal. When the 
granules "do not touch", each moves independently of the other (i.e., their potential of interaction is zero beyond a 
certain distance between their centers). When they do "push one into the other", there is a repulsive linear force 
between their centers with force constant k; the force vanishes once the particles separate again. One or both of the 
granules are subject to friction, which may be kinetic (i.e. constant, see Sec.|nJ or hydrodynamic (i.e. proportional 
to velocity, see Sec. IIII|) . In order to initiate energy transfer between the granules, one can, for example, introduce 
some precompression in the system, or one can impart one of the granules an initial velocity toward the other. We do 
the latter: initially the granules are assumed to be just touching, the first granule has initial velocity vo toward the 
second, and the second is at rest. We are interested in the velocity of the two granules and the total energy in the 
system immediately after the collision. Note that in the absence of friction the collision is elastic so that the energy 
and momentum are conserved. In this case the first granule simply ends up at rest and the second with velocity Do, 
i.e., there is no backscattering. Of course there is backscattering, even in the absence of friction, in a longer chain. 

Even these simple systems present a broad range of outcomes, some more interesting than others, and not all of 
which will be covered in detail. For example, one or the other or both granules may stop moving altogether before 
the collision is over. However, we are particularly interested in situations where the second granule is still moving 
forward after the collision (propagation) and the first granule is moving backward (backscattering) . Our analysis will 
focus on these outcomes and the conditions that lead to them. 

In Sec.[n]we present and solve the two-particle model with kinetic friction. Section llTll presents a similar calculation 
for the model with hydrodynamic friction. A summary of results, and a brief outlook of future work, are gathered in 
Sec.HVI 

II. TWO PARTICLES WITH KINETIC FRICTION 

In this section we examine granule collision when one or both of the granules are subject to kinetic friction with 
friction coefficient fi. Kinetic friction arises from contact forces when a solid body moves along a solid surface, and 
depends on the velocity only in that its direction is always opposed to the velocity. Two cases are presented analytically 
(at least in part). They turn out to exhibit rather different behaviors, which we anticipate here as a point of reference. 
In one, only the second granule is subject to friction. We find that for a certain range of friction coefficient values there 
is both propagation and backscattering, that the velocities of the forward moving second granule and the backward 
moving first granule are monotonic functions of the friction, but that the total energy dissipation during the collision 
is nonmonotonic. In the second model only the first granule is subject to friction. Here we again find a range of 
values of the friction coefficient for which there is both propagation and backscattering. Now, however, the granule 
velocities are nonmonotonic functions of the friction while the energy loss is a monotonic function. Results for the 
more general case of both granules subject to (possibly different) friction are presented numerically. In this way we 
cover the entire range of "frictional disorder", as might occur if the granules are made of different materials. While 
most of our discussion revolves around the case n = 2, some results for n — 5/2 are presented. Also, as a backdrop we 
keep in mind the most extreme case of a hard interaction, a hard sphere (n — + oo). The interaction is instantanous 
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in this case and leads to elastic collisions in which energy and momentum are conserved. The first granule stops and 
the second granule moves forward. This comment emphasizes the observation that only collisions of finite duration 
can lead to a combination of backscattering and propagation in two granules. 

We rescale variables (see Appendix [SJ so that we can set the coefficient a = 1 in Eq. regardless of the value of 
n. This rescaling allows presentation of results in terms of a single parameter, the scaled friction. 



A. Frictionless First Granule 



Consider first the case of a frictionless first granule. The movement of the second granule will then only start if the 
clastic compression exceeds the frictional force, > 77, and so at first only the first granule moves and its motion 

is governed by the equation of motion 

x 1 (t) = -[x 1 (t)] n -\ (4) 

The initial condition for this problem is xi(0) = and xi(0) = 1. The solution for n — 2 is easily found to be 

x 1 (t)=sm{t), (5) 

and is valid until the time to at which xi(to) = 77, that is, until to = arcsin(?7). 

At time to the second granule starts moving as well, and the system is now governed by the two coupled equations 
of motion 

Xl (t') = -[x 1 {t')-x 2 {t 1 )] (6a) 



x 2 (t') 



-r?- [x 2 (t') - xi(t')} 



(6b) 



with the initial conditions: x\(t' — 0) = 7/, ±i(t' = 0) = y/1 — rj 2 , and ^(t' = Q) = X2(t' = 0) = 0. Here we have 
defined t' = t — tg- The solution is found to be: 



-r?(-3 + t' 2 ) + 2(1 - ri 2 ) 1/2 t' + r)cos(V2t') + \/2(l - r; 2 ) 1/2 sin(v / 2t') 



(7a) 



x 2 (t') 



T)(-l + t' ) - 2(1 - n 2 ) 1/2 t' + 77C0s(V2i') + 72(1 - ry 2 ) 1/2 sin(v / 2t') 



(7b) 



If 77 is too large, the second granule will stop before the two granules lose contact. For sufficiently small 77 this 
does not happen and the above solutions then remain valid until contact is lost. The time t' = t\ when this occurs is 
the time at which Ax(ii) = x\(t\) — X2(t\) vanishes, that is, when t\ = The velocities of the granules at the 

moment of separation t = to + t\ are 



77 7T 
2V2 



» 



((l-ry 2 ) 1/2 



2V2J 



(8a) 
(8b) 

The values of 77 for which this description is valid are those for which t/j > 0, that is, for < 77 < (1 + 7r 2 /8) -1 / 2 = 
0.669 . . . Beyond time t = to + 1\ each granule continues on its course, the first forever and the second until friction 
brings it to rest. Backscattering and propagation thus occur when 77 < 0.669 ... As a function of 77 in this range, the 
magnitude of the backward velocity increases monotonically while that of the forward velocity decreases monotonically. 
The backscattering is shown in Fig. ^ which we exhibit in order to stress the similarity in the behavior of the n = 2 



,(«) 



and V2 are 



potential and that of the Hertz potential for spheres. As n increases the range of 77 for which both 
nonzero shrinks, until it shrinks away completely for a hard sphere potential. 

Perhaps the most dramatic outcome of this calculation is the fact that the friction on the second granule causes the 
first granule to move backward (recall that if both granules are frictionless the first one simply stops). In other words, 
the friction on the second granule is responsible for backscattering. This behavior, strange at first glance, occurs 
because the first granule "stops too soon" , that is, it stops before the spring has had time to decompress completely. 
Indeed, if the friction on the second granule is so strong that it does not move at all, the first will move backward 
with velocity v\ — —1. Our model lies between the limiting cases of zero and very large 77 (-1 < v[ s) < 0). 



FIG. 1: Velocity of the first (frictionless) granule at the end of the collision as a function of friction. Solid curve: n — 2.5 (Hertz 
potential for spheres). Dashed curve: n = 2. 




FIG. 2: Ratio of the total energy at the moment of separation to the initial energy as a function of friction. Solid curve: n = 2.5 
(Hertz potential for spheres). Dashed curve: n = 2. 



In Fig. |3 we show the ratio of the total kinetic energy of the system at the end of the collision to the initial energy, 



E 



1 v(i-v 2 ) 1/2 - f* 3 1 

V2 V 4 



(9) 



In the rescaled variables, E = 1/2. Curiously, there is a region where this ratio increases with increasing friction 
coefficient. The reason for this behavior is that when friction increases beyond a certain point (beyond the minimum 
in the curve), relatively more energy is transfered to the first granule (as opposed to being simply dissipated) by the 
end of the collision. Again we also show the numerical results obtained for the Hertz potential with n = 5/2 (solid 
curve) and comment on the similarity with the n — 2 curve. 



B. Frictionless Second Granule 



Next we consider the case of a frictionless second granule. The first granule now moves forward, eventually stops, 
and perhaps moves backward, while the second granule moves forward. 

As before, the granules are initially just touching, the first granule having an initial velocity vq = 1. While the first 
granule is moving forward the system is governed by the equations 



Xl{t) 



-1] - [xi(t) - x 2 {t)} 



(10a) 



x 2 {t) = - [x 2 (t) - XI (t)] 



(10b) 



FIG. 3: Comparison of the small-77 analytical solution (dotted) and the numerical solution (solid) of the stopping condition 
vi(to) = 0. The dashed line is the time it would take the granules to separate if the solutions lllat and All bi t remained valid 
until then. The fact that it lies above the other curves indicates that the first granule stops before the collision is over. 



Explicit integration yields 



Xl(t) 



x 2 {t) 



1 r 

4 

1 

4 



-77 + 2t - r/t 2 + rj cos(V2t) + V2 sin(v / 2t) 
77 + 2t - rjt 2 - r\ cos(\/2i) - y/2 sin(v^t) 



These equations hold as long as the velocity v\(t) = Xi(t) > 0. This velocity is given by 



ui(t) 



2 - 2-qt + 2 cos( V2t) - V2rj sin( V2i) 



(11a) 



(lib) 



(12) 



The time to at which v\(to) = cannot be found explicitly in general, but the problem may still be treated analytically 
for small or for large 77, or numerically for arbitrary r). 

If 77 1, we can expand «i(t) around t = (the solution for rj = 0). We find that 



tr. 



v/2 



V2 



1/2 



+ 0(t? 3/2 ) 



(13) 



In Fig. we compare this analytical approximation with the numerical solution. The small- r\ approximation is clearly 
very good for a rather large range of r\. Therefore, we proceed using this approximation. In the figure we also show 
the time that it would take the collision to be over if the solutions (|llafl and ljllb|) would remain valid, i.e., the 
solution t' of the equation x\{t' ) — X2(t' Q ). This time, shown as a dashed line, is always above the time at which the 
first granule stops, indicating that the first granule necessarily stops before the collision is over. This persists even 
for large r\: as 7/ — > 00 the stopping time of the first granule goes to zero as I/77, while the time for the collision to be 
over if the above solutions would remain valid goes to zero as 2/77. 

When the first granule stops, the equations of motion are no longer Eqs. (|10af) and l|10b|) . The initial conditions for 
the new set of equations are the velocities and positions found above at time t — tg. Within the small- 77 approximation 
we use 



Xl(t 



2V2ir - (tt 2 + 4)77 



(14a) 



X2(t ) 



2V21T - 4(2 3/4 )y/¥fj - (tt 2 - 4)77 



(14b) 



Vi(to 



(14c) 



V2(to) 



V2TTT] 



(14d) 



as the initial conditions. 
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FIG. 4: Comparison of time required for the two granules to separate calculated numerically (solid) and analytically for small 
rj (dotted). The time required for the first granule to stop according to the equations of motion 116aH and dl6bl) is shown as a 
dashed curve. Therefore, if rj < 0.22 the first granule is still moving backwards at the end of the collision. 



The equations of motion at this point depend on whether the first granule simply remains stationary or starts 
moving backward. The latter occurs if the elastic force is greater than the friction, that is, if 



[xi(to) - x 2 (t )} 



-v-v 



cos^ 1 / 4 ^) + V2 sin(2 1 / 4 07?) 



> V- 



(15) 



Therefore, for 77 smaller than a critical value (77 = 0.439 . . .) the first granule will start moving backward, and we 
confine our analysis to this case. The equations of motion then are (note that since the first granule is moving 
backward, the friction is also reversed): 



Xi 



(f) = v -[x l (t')-x 2 (t')\ 



(16a) 



£a(0 = -[a:a(0-si(0] ( 16b ) 

where t' = t — to- The initial conditions are those given in Eqs. JUaJ-JIUJ. 

The solutions of this set as a series in 77 can easily be exhibited explicitly but are not very instructive. Two cases 
must be distinguished: 1) The first granule is still moving backward when the collision ends; and 2) The first granule 
stops before the collision ends, at which point the equation of motion (|16a|) no longer holds. In either case, from 
the series solutions one can calculate the time t' = t\ at which the collision ends as the solution of the equation 
x\{t') = X2{t'). If we assume that the first granule is still moving we find 

tl =(^ 1/ °-T, + Otf'*). (17) 

This solution is shown as the dotted curve in Fig. 0] as is the numerically obtained exact solution (solid line) . The 
approximation i|17|) begins to seriously deviate from the exact solution at around 77 ^ 0.1. The dashed curve indicates 
the time at which the first granule stops before the end of the collision if Eq. i|16a|) remains valid until it does. This 
analysis is therefore appropriate and leads to backscattering (and propagation) if 77 < 0.22. 

In Fig. we see that as 77 increases the final backward velocity of the first granule is nonmonotonic, that it, is 
increases from zero to approximately 3% of the initial velocity, and then decreases again. In this model, we thus 
observe that there is backscattering for small 77. For 77 = and for 77 greater than a certain value (> 0.22), the first 
granule is at rest when the collision ends and there is no backscattering. Therefore, the backscattering in this model 
happens because the first granule (as opposed to the second granule in the previous example) "stops too soon" (i.e. 
before the collision is over), but not "too late" (i.e. when the compression is still strong enough to overcome friction). 
The required balance leads to the nonmonotonic behavior seen in Fig. Even though the backscattering is small 
(e.g. compared to that in the first model), it is conceptually striking. However, the total energy of the system at the 
end of the collision is a monotonically decreasing function of 77, in contrast to the behavior in the previous case shown 
in Fig. H 



FIG. 5: Final velocity of the first granule. 



C. Granules with Arbitrary Friction 



While the results presented above deal with the extremes of frictional disorder, it is helpful to present a broad 
brush of graphical results for arbitrary combinations of 771 and 772, the friction parameters for the two granules. This 
is also an opportunity for re-assertion that these broad brush results are essentially the same for the n — 2 and the 
n = 5/2 Hertz potentials for two granules. Since the graphical results are necessarily numerical, we present them for 
the n = 5/2 potential. They are qualitatively indistinguishable from those for n = 2. 

In Fig. we show in gray scale the magnitudes of v[ s ^ (first panel), (second panel) and energy ratio E^/Eq 
at the end of the collision in the (771,772) regime where there is simultaneous backscattering and propagation. Darker 
shading represents a higher speed (first and second panels) and a higher energy (third panel). Figure ^ is associated 
with an upward trajectory in the first panel at 771 = 0. The monotonically increasing backscattering seen in Fig. 
corresponds to the upward darkening. Figure [21 is associated with the same upward trajectory, now in the rightmost 
panel. The nonmonotonicity is seen in the lightening followed by darkening. The nonmonotonicity of Fig. [5] is 
associated with the shading in the first panel associated with the horizontal trajectory along 772 = 0, captured in the 

curvature of the iso- velocity lines. The monotonicty of discussed in the (771,0) and (0,772) cases is evident in the 
middle panel. 

This rather complex behavior arises from the interplay of energy dissipation and energy transfer during a collision. 
Thus, for example, increasing friction on the first granule causes it to stop sooner. If it stops too soon ("high 771"), 
the spring is not sufficiently compressed to overcome friction and it may happen that the whole process simply stops 
then and either one or the other or both granules do not move. On the other hand, if it stops too late ("low 771") then 
the energy transfer to the second granule may cause it to move forward, but there may not be sufficient energy to 
cause backscattering. Or backscattering may occur but not forward motion. Between these situations ("intermediate 
771"), both granules may be in motion at the end of the collision. 



When dissipation is of hydrodynamic origin, the dissipative force is proportional to the particle velocity, — jx. 
Hydrodynamic friction arises from the motion of a body against a lubricated surface or in a liquid or gaseous medium. 
We find that the cases of equal and slightly unequal friction coefficients can be solved analytically and cover essentially 
the whole spectrum of possible behaviors. 

The scaled equations of motion (see Appendix^) for the potential with n = 2 are 



III. TWO PARTICLES WITH HYDRODYNAMIC FRICTION 



xi(t) 



7i±i - [xi(t) - x 2 (t)] 



(18a) 



3*2 (*) 



72± 2 - [aiaW - xx{t)] 



(18b) 



Even though these equations can be integrated exactly, the time required for the two granules to separate can, in 
general, only be calculated numerically. Here we consider the case 71 f=a 72, for which we can find an approximate 
analytic expression for this time. Thus, we set 7$ = 7(1 + Si) with 81 = —62 = S and solve the problem for small \5\. 
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FIG. 6: Backscattering velocity — (first panel), forward velocity (second panel), and energy ratio /Eo (third panel) 
at the end of a collision for the Hertz potential for spheres. In the first two the gray scale shows darker coloring for higher 
speeds. In the third panel darker coloring indicates higher energy. The curves in the first two panels are iso-velocity lines, and 
in the third panel iso-energy lines. 

Defining x± = x\ ± X2, we can obtain from Eqs. I|18|) the associated Laplace transforms X±(s): 

5jsX_ (s) + (js + s 2 ) X + (s) = 1 (19a) 

(2 + 7 s + s 2 ) X_(s) + 5jsX + (s) = 1. (19b) 
The solution of these equations to first order in 8 is 

1 5-i 

X+(s) = . rn ' . -r (20a) 

1 5i 

X -{s) = 7, 7 s-7 -7- — - — 7 rr. (20b) 

w 2 + s (7 + a) 7 + a) 2 + s (7 + s v ; 



FIG. 7: Final velocity of the first granule when 71 > 72, S — 0.1 (dashed), 71 = 72, 5 — (dotted), and 71 < 72, 5 = —0.1 
(solid). 



The time functions x±(t) are the inverse Laplace transforms of X±(s). Furthermore, the time to at which the granules 
separate can be calculated by solving the equation X-(t = to) = 0. To zeroth order in S (that is, for 8 = 0), we have 

x_(t ) = e < to/2 — — = 0, (21) 



where uj = y8 — 7 2 /2. If 7 2 > 8 the collision never ends. If 7 2 < 8 then it ends at time to — ir/ui. To calculate 
the separation time to first order in 5, we write t% = to + e/w, expand X—(t%) = to first order in e, and solve the 
resulting equation for e, to obtain 

e =-^(l + e -WV). (22) 
The final velocities of the granules at the end of the collision to first order in 5 are: 

v[ s) = - (l - e (7r/2uJ ^j e ^ 77r/w + (l + 2e 77r/2w + e 77r/ ") e " 377r/2 " (23a) 



4 S) = i (l + e^' 2 ") e-^^ + ^ [Slu 2 + (Uu; 2 - 7 2 + 8) e 7 ^ 2w + (4co 2 - 7 2 + 8) e^^] e ~^/^. (23b) 



We immediately see that the backscattering is greater if 72 > 71 (S > 0) than if 71 is the larger friction coefficient. 
This behavior is similar to that which occurs when we have kinetic friction. Thus, again a force in the second granule 
is more effective in producing the backscattering (as expected) . In Fig. we present the velocity of the first granule 
at the end of the collision as a function of 7 for different values of 6. Once again we find (for each value of 8) an 
optimal scaled dissipation constant for which the backscattering is greatest. We note that the total energy at the end 
of the collision is a monotonically decreasing function of the average dissipation parameter 7. 



IV. SUMMARY AND OUTLOOK 



We have presented an overview of the dynamics of two granules subject to repulsive but not restoring forces (as 
in granular materials), and to kinetic or hydrodynamic friction. In particular, we explored parameter regimes that 
would give rise to backscattering and propagation at the end of a collision, with a view toward future studies of 
chains of granules subject to friction. We explored the consequences of "frictional disorder" in the two-granule system 
and found dramatically asymmetric behavior, that is, backscattering and propagation are affected differently by the 
friction on the first and second granules. For example, while friction that is "too high" on either granule will lead to 
situations where backscattering or propagation or both do not occur, the range of friction parameters that do lead 
to the occurrence of backscattering and propagation is considerably greater for the second granule than for the first. 
While this is perhaps not surprising, the nonmonotonicities in the final velocity and energy as a function of friction 
point to the interesting dynamical asymmetries in this problem. While most of our results have been calculated for a 
potential with n — 2 [cf. Eq. l(T])]. we have argued and shown numerically that the two-granule results are qualitatively 
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similar for other values of n (albeit not quantitatively equal). They therefore provide some insights into the dynamics 
involving more realistic potentials. 

Finally, we stress that a chain of many granules differs in important ways from the two-granule system. For example, 
with n > 2 a frictionlcss chain may support solitons while the n — 2 chain does not [3J. Also, due to the interactions 
with further granules, there may be backscattering in a chain even for parameters that do not lead to backscatteringin 
the two-granule system. On the other hand, since pulses are known to remain very narrow in granular chains |2Lblul> 
the analysis of the dynamics of two granules provides important insights into the behavior of longer chains and other 
larger arrays [l8j . 
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APPENDIX A: SCALING 



The equations of motion with kinetic friction are in general of the form 

mj/i(r) = ±mg - a[y 2 (r) - ^(r)]™" 1 (Al) 

where the t/'s are displacements, r is the time, n is the exponent in Eq. lJT|l. fJ-i is the friction coefficient, g is the 
gravitational constant, and i,j = 1,2. If particle i is frictionless, the first term on the right is absent. The initial 
conditions are t/i(0) = and j/i(0) = or vq. Wc define new variables Xi and t via the relations, 

yi = Ax i: r = Bt, (A2) 

in terms of which the equations of motion are 

mA . 
~W 

The choices 



— xi(t) = ± m - oA*- 1 - ^(t)]"- 1 . (A3) 



/ 2 \ V n i / 2 \ !/" / 2 \ V™ 

A= Ml) , B=l(^) , „= J* (ML) (A4) 
\ a J v \ a J towq \ a J 

lead to the scaled equations of motion used in the text, e.g. 

Xi(t) = ±m - \xi(i) - x 3 {t)] n '\ (A5) 

with the initial conditions 2^(0) = and ii(0) = or 1. 

With hydrodynamic friction the equations of motion are of the form 

my\{T) = -7^(t) - a[ yi { T ) - y^r)]"- 1 (A6) 



and, in terms of the new variables, 



A . _ A 

Again we choose A and B as above, and set 



Xi(t) = -7ii*i(*) " aA n - x [ Xi {t) - x j (t)] n ~ 1 - ( A7 ) 



/ 2 \ 1/™ 

7* = ^ (A8) 
itivq \ a J 

to obtain 

*,•(*) = - [Xi(t) - x^t)} 11 - 1 (A9) 

with initial conditions as above. 
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